Subspaces that Minimize the Condition Number of a Matrix
نویسندگان
چکیده
We define the condition number of a nonsingular matrix on a subspace, and consider the problem of finding a subspace of given dimension that minimizes the condition number of a given matrix. We give a general solution to this problem, and show in particular that when the given dimension is less than half the dimension of the matrix, a subspace can be found on which the condition number of the matrix is one. 1 The problem Suppose A ∈ R and V ⊆ R is a subspace with dimV = k ≥ 1. We define the maximum gain (minimum gain) of A on V , as Gmax = sup x∈V, x 6=0 ‖Ax‖ ‖x‖ , Gmin = inf x∈V, x 6=0 ‖Ax‖ ‖x‖ , respectively, where ‖ ‖ denotes the Euclidean norm. When A is nonsingular, we define its condition number on the subspace V as κV(A) = Gmax/Gmin. The condition number of A on any one-dimensional subspace is 1, and its condition number on V = R is the (usual) condition number of A, which we denote κ(A). The condition number on any subspace is between 1 and κ(A). If κV(A) = 1, we say that A is isotropic on V , since its gain ‖Ax‖/‖x‖ is the same for any nonzero vector x ∈ V. In this note we address the following problem: Given a nonsingular matrix A ∈ R, and k ∈ {1, . . . , n}, find a subspace V ⊆ R of dimension k which minimizes κV(A). The number κV(A) is a measure of the anisotropy of the linear function induced by A, restricted to the subspace V , so our problem is to find a subspace of dimension k on which A is maximally isotropic.
منابع مشابه
An open letter concerning Subspaces that Minimize the Condition Number of a Matrix
We define the condition number of a nonsingular matrix on a subspace, and consider the problem of finding a subspace of given dimension that minimizes the condition number of a given matrix. We give a general solution to this problem, and show in particular that when the given dimension is less than half the dimension of the matrix, a subspace can be found on which the condition number of the m...
متن کاملStructured Condition Numbers for Invariant Subspaces
Invariant subspaces of structured matrices are sometimes better conditioned with respect to structured perturbations than with respect to general perturbations. Sometimes they are not. This paper proposes an appropriate condition number cS, for invariant subspaces subject to structured perturbations. Several examples compare cS with the unstructured condition number. The examples include block ...
متن کاملPerturbation Bounds for Isotropic Invariant Subspaces of Skew-Hamiltonian Matrices
Abstract. We investigate the behavior of isotropic invariant subspaces of skew-Hamiltonian matrices under structured perturbations. It is shown that finding a nearby subspace is equivalent to solving a certain quadratic matrix equation. This connection is used to derive meaningful error bounds and condition numbers that can be used to judge the quality of invariant subspaces computed by strongl...
متن کاملPreconditioned Generalized Minimal Residual Method for Solving Fractional Advection-Diffusion Equation
Introduction Fractional differential equations (FDEs) have attracted much attention and have been widely used in the fields of finance, physics, image processing, and biology, etc. It is not always possible to find an analytical solution for such equations. The approximate solution or numerical scheme may be a good approach, particularly, the schemes in numerical linear algebra for solving ...
متن کاملA Framework for Adapting Population-Based and Heuristic Algorithms for Dynamic Optimization Problems
In this paper, a general framework was presented to boost heuristic optimization algorithms based on swarm intelligence from static to dynamic environments. Regarding the problems of dynamic optimization as opposed to static environments, evaluation function or constraints change in the time and hence place of optimization. The subject matter of the framework is based on the variability of the ...
متن کامل